Coherent transport in multi-branch quantum circuits 
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An open multi-branch quantum circuit is considered from the viewpoint of coherent electron or 
wave transport, both with and without intrinsic disorder. Starting with the closed system, we give 
analytical conditions for the appearance of two isolated localized states out of the energy band. In 
the open system, using the method of the effective non-Hermitian Hamiltonian, we study signal 
transmission through such a circuit, with an important result of a long lifetime of localized states. 
When the average level width becomes comparable to the mean level spacing, the super-radiant 
transition occurs. In the case of on-site disorder we find an analytical estimate, confirmed by 
numerical data, for the robustness of the isolated states and their role in transport processes. 
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INTRODUCTION 

The development of quantum informatics requires bet- 
ter understanding of the general problem of quantum sig- 
nal transmission through discrete structures of interact- 
ing quantum elements, such as quantum dots or 
molecular 0, Q and Josephson junctions [f|. There is 
growing theoretical and experimental interest in arrange- 
ments more complicated than a simple one-dimensional 
chain, including Y- and T-shaped structures M,l8|, tetra- 
hedral qubits [9j, connected benzene rings [lOj, crossed 
chains [H[ , two- and three-dimensional lattices [HI, EH , 
as well as graphs with the violation of time-reversal in- 
variance at the junctions, supposedly with the aid of mag- 



netic fields [14j . In all cases, one has to deal with a quan- 



tum system with intrinsic stationary states that become 
unstable when the system "opens" to the external world 
as a part of a transmission network. 

The method of the effective non-Hermitian Hamilto- 
nian, borrowed from nuclear physics 15], is a powerful 
and efficient tool for theoretical analysis of open quantum 
systems, as was shown by applications to one-dimensional 
structures covering both regular and chaotic 

internal dynamics; see the recent review of the method 
and its various adaptations 



21] 



In this work we con- 
sider the M-branch circuit in the form of M > 2 one- 
dimensional tight-binding chains with a common vertex 
at the central point. This system is a simple discrete ex- 
ample of quantum mechanics in a non-trivial space with 
self-crossings, singular points or surfaces where quasi- 
bound states (evanescent waves) may emerge even for 
open boundary conditions. Our interest is not only in 
the structure of the energy spectrum and eigenfunctions 



for closed samples, which was earlier discussed for differ- 
ent purposes, see for example, the applications to Bose- 
Einstein condensation and temperature effects in optical 
networks 22j, but mainly in the transport characteris- 
tics for the case when the chains are connected to an 
environment. 
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Figure 1: (Color online) The M-branch circuit with the coupling 
at the origin. 



MODEL 

We start with the discussion of the closed system and 
its energy spectrum. In our model (see Fig. [T]), each 
of the M branches consists of N a sites (a — 1,...,M) 
along which a particle/wave can propagate through the 
structure. The hopping amplitude v a between nearest 
neighbor sites is constant within each branch, while the 
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coupling to the origin is given by t a v a . Thus, the Hamil- 
tonian of the closed system is 



M Na — l 



H=^ v a (\a,n)(a,n+l\ + \a,n + l)(a,n\ 



a—l n—1 

M 
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In what follows, unless stated otherwize, we set v a = v 
and N a = N, so that there are K = MN + 1 sites in 
total, including the vertex. The Schrodinger equation 
for a stationary state with energy E reduces to a set 
of algebraic equations for the site amplitudes C° , where 
a = 1, . . . , M, and n labels the sites in each chain, 

A I 

Y^tavCt = EC : vC%+t a vC = EC?; 

a=l 

«(C£-i + C a n+l ) = EC a n for 2 < n < N , (2) 
with the boundary conditions C^ +1 = 0. 

SPECTRUM AND EIGENFUCTIONS 

For Co 7^ 0, the a-dependence of the amplitudes is 
given by C° = vt a X n with n = 1, . . . , N, while Co and 
the a-indcpcndcnt amplitudes X n obey a tri-diagonal ho- 
mogeneous set of N + 1 linear equations. The secular 
equation for E ^ ±2v has the form, 

(E 2 + Ee+ - v 2 Q)e N - (E 2 + Ee_ - v 2 Q)e^ = , (3) 

with the control parameter Q , 
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M 



(4) 



which generalizes the results found in Refs. 23 1. 

In order to find the structure of the corresponding 
eigenfunctions from Eq. @, we obtain, after lengthy cal- 
culations, the following relation, valid for any finite N: 



Ga 
n+1 



Cff 



— Bn(Q,N), n=l,...,N-l, (6) 



where we have defined the boundary factor, 



B n (Q,N) 



i-(Q-i) 



-(N-n) 



1-(Q-1 



-(N+l-n) 



(7) 



which differs from 1 only for n ~ N. Neglecting for the 
time being this boundary effect (which however will be 
crucial for an open model, see below), we obtain 



Co 



Q-i 



-(n-l)/« 



n = 1, 



Here ^ is the localization length, 



1 



HQ - 1) 



■ ,N; Q>2. (8) 



(9) 



of the two states outside the energy band, with the peak 
located at the origin of the circuit. Note that these states 
are spread over all branches in proportion to the cou- 
plings t\. 

For Co = 0, one has standard extended Bloch states 
in each branch. However, the M solutions are linked 
through the first equation of ((2]), yielding M — 1 inde- 
pendent degenerate states for each Co = eigenvalue. As 
a result, for the symmetric case of equal branch lengths 
N a = N, we have two isolated localized states, sets 
of M — 1 extended degenerate Bloch states with Co = 0, 
and (N—1) extended non-degenerate Bloch states with 
C ^ 0, altogether K = NM + 1 states. 



and e±(E) = (-E±\ / E 2 - 4v 2 )/2. Note that for even N 
the value E = is also a solution of Eq. ^ . An analysis 
of Eq. shows that the energy spectrum consists of a set 
of MN — 1 eigenvalues within the energy band \E\ < 2v 
(with eigenfunctions extended over the M branches of the 
circuit), and two additional eigenvalues with \E\ > 2v, for 
sufficiently large M or sufficiently large couplings t a . 

We first consider the two states with the eigenvalues 
outside of the energy band (sometimes such states are 
called "hidden"). These states turn out to be strongly 
localized at the origin of the circuit provided the param- 
eter Q is large enough. In the limit N ^> 1, the energies 
of these states can be found from Eq. ([3]) by considering 
its largest and smallest roots. For Q > 2 one obtains 



floe I — 



Q 



(5) 



COUPLING TO CONTINUUM 

Coming to the main goal of our study, now we consider 
the same circuit coupled to the environment by attaching 
the last site of each branch to an external channel, simi- 
larly to what has been done in Refs. [H [H, 0, 0] . The 
open system is described by an effective non-Hermitian 
Hamiltonian [H|,[2l|, 

M 

^ = H -1 7 W; W = Y,A a N \a,N)(a,N\A%*. (10) 

a=l 

Here H is given by Eq. (p}, while the coupling to the con- 
tinuum via the ends of branches is characterized by the 
parameter 7 and the matrix W. The matrix W is con- 
structed out of the transition amplitudes between 
intrinsic states |a, N) and continuum states \a, E). We 
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set = 1, so that the strength of the coupling is con- 
trolled entirely by the parameter 7. 

The lifetime t\ oc of the two localized states can be 
estimated via the imaginary part — ri oc /2 of the corre- 
sponding eigenvalues £\ oc of the complex Hamiltonian, 
t\ oc = rr c . For small 7, the resonance width is deter- 
mined by the spatial overlap of the localized states with 
the edges, 



M 

7^|(a,iV a |^oc)| 2 

a=l 



(11) 
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Taking into account Eqs. ([HI El El) with 
for N, Q > 1 we 



see Refs 

I Co 1 2 found from the normalization 
have 



"Hoc = 2 



(Q-i 



.N+l 



7<3(Q-2) 
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7 



N-l 



(12) 



Considering now the extended states inside the energy 
band, two different regimes can be distinguished [lEl [l6[ 
as a function of 7. At weak coupling, all these states are 
similarly affected by the continuum coupling and acquire 
widths proportional to 7. For large 7, only M "super- 
radiant" states have a width proportional to 7, while 
the widths of the remaining ("trapped") states fall off as 
I/7. In order to find the critical value of the parameter 
7 corresponding to the super-radiant transition, we ana- 
lyze the average value (r) of the (MN + 1) — M narrow 
widths as a function of the rescaled coupling 7/11. At a 
critical value 7 cr , the average width (r) peaks and begins 
to decrease. 
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Figure 2: (Color online) The dimensionless average width (T) as a 
function of J /v. We consider a) M = 4 and N = 70 with t a = 1 for 
all branches (red triangles); b) M = 4 and N = 70 with different 
couplings t\ = 200, t2 = 100, t$ = 20, and ti = 50 (blue squares); 
and c) M = 4 with different numbers of sites, N\ = 55, = 
80, N 3 = 58, N4 = 87, and all t a = 1 (green circles). The dashed 
line corresponds to the average over all widths, while the symbols 
are obtained by averaging over the K — 4 smallest widths, where 
K = 281 is the total number of sites in all cases. 

One can evaluate 7 cr using the following criterion 
lil 20j j : the transition occurs when (r) becomes of the 



for the closed system. This is particularly easy for the 
special case of equal coupling t a — t and equal num- 
ber of sites in each branch N a — N; in this case it is 
convenient to define an effective mean level spacing as 
D c g w 4v/(2N + 1) ~ 2v/N, which takes into account 
the level degeneracy. On the other hand, the average 
width is given by (r) M~//(MN + 1) ~ 7/A, so that 
7 cr = 2v independently of M, N, and t a . This result is 
numerically confirmed in Fig. [2j where it is shown to be 
valid even in the more general case of different coupling 
t„ and different number of sites N„ . 



TRANSMISSION 

Maximum transmission occurs at the super-radiant 
transition, as happens in one-dimensional chains [20j |. 
Moreover the above analysis allows one to understand 
generic properties of the transmission between different 
branches. The data in Fig. [3] demonstrate the energy de- 
pendence of the transmission coefficient T a b between two 
channels b and a. One can see that with an increase of 
7 the resonances corresponding to the energies from the 
bulk of spectrum begin to overlap, in contrast to the two 
very narrow resonances located out of the energy band. 
These quasi-bound resonances remain extremely stable 
even for a very strong coupling to continuum. 
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Figure 3: (Color online) Transmission coefficients T12 as a function 
of energy, for M = 4, ti = 3, ti = 2, and t 3 = t 4 = 1. Full red 
curves stand for 7 = 0.6, while dashed blue curves correspond to 
7 = 5. a) symmetric case N = 4; b) asymmetric case iVi = 4, ^2 
■ JV3 = 3, -/V4 = 5. Horizontal dotted lines are the theoretical 
results given by Eq. (1 1 5 P and its generalization to unequal branch 
lengths. 

Of special interest is the value of the transmission co- 
efficient T a b for the resonances outside the band. Since 
these resonances are very narrow, one can estimate 
T a b(£ioc) [llllill in terms of the edge components of the 
localized state tf>i oc in channels a and 6, 



order of the mean level spacing D of the Hamiltonian 



(a, N a \^ oc )^f(^ loc \b, N b )^f 



(i/2)r loc 



(13) 
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which can be computed from Eqs ((6l [7J : 



N a -1 



\{a,N a \<k oc )\ 2 = \C \ 2 t 2 a e- N ^ [] B n (Q,N a ). (14) 



In the special case of equal-length branches this result 
assumes the particularly simple Hauser-Feshbach form: 



Tab(£\oc) 



At 2 ti 



(15) 



This relation is in good agreement with the data of 
Fig. [3^, as indicated by the dotted horizontal line. The 
value of T a b{£\ oc ) for narrow resonances depends on the 
hopping elements t a only, and is independent of the cou- 
pling strength 7 for equal values of the couplings, 7a = 7- 
The numerical data also indicate that, for equal- length 
branches, the maximal value of the transmission coeffi- 
cient for \E\ < 2 is given by the same expression (TT51) . 
while this does not happen when branches have differ- 
ent length, see Fig. [3]b). As 7 — > 0, the lifetime of the 
localized states becomes very long. 



INTRODUCING DISORDER 

A key consideration in important practical applications 
is the influence of disorder. For this reason we study the 
circuit with diagonal disorder by adding the term V, 



V = ^2 2J e a,n\a, n) (a, n\ 



(16) 



to the non-Hermitian Hamiltonian (fTU)) . Here the site 
energies e a ,n are random numbers uniformly distributed 
in the interval [-W/2, W/2]. According to the theory of 
disordered systems, for N — > 00 all states within the en- 
ergy band become exponentially localized with a localiza- 
tion length cx (v/W) 2 . For the isolated localized states, 
a weaker dependence on disorder can be expected, and 
it is more convenient to define a localization length i\ oc 
through the inverse participation number, 



critical value can be estimated by assuming that it cor- 
responds to the intersection between the gap of size A 
emerging due to disorder around the localized states and 
the bulk of the spectrum of width 2(2v + W/2). One can 
estimate the value of A from the relation, 



A 



{lp\oc\V\ll)\ oc 



12 



7-1 
loc ' 



(19) 



where the average is taken over sites n and over the dis- 
order. The critical value W cr for large N is thus ob- 
tained by equating the half- width of the density of states 
(2v + W/2) to the minimal possible energy of the isolated 
states due to random fluctuations: 



2v 



W 
~2 



(20) 



which yields 



Wcr =2V 



Q 



- 2 



\ 



\ 



(Q-2) Q(Q-2) + J2^ 1 t 



12Q(Q-iy 



J 



(21) 

Despite the cumbersome appearance of Eq. (|2"Tj) . it ad- 
mits two interesting limits: assuming for simplicity t a = t 
for all a, it is easy to see that for Q w 2, W CI /v ~ 
(Q — 2) 2 /2, while for Q — > 00, one has W cr /v ~ ky/Q, 
where the constant k = 2/(1 + y/(M + 1)/(12M)) de- 
pends only on the number of branches, and k w 1.55 
for large M. Note that it is impossible to have delocal- 
ized isolated states along with localized states inside the 
band for any disorder W. The expression ()21j) is com- 
pared with numerical data for W a defined as the point 
where the localization length £\ oc of the isolated localized 
state begins to decrease with an increase of disorder. The 
data in Fig. [4] show quite good agreement with this es- 
timate over many orders of magnitude with no fitting 
parameters. 



M N 



l(o|Vioc}| 4 + ]TEK a '^ioc) 



(17) 



a— 1 n— 1 



Without disorder and for large N, the value of l\ oc for 
the isolated states can be estimated as 



t4oc 



4Q(Q-1) : 



(Q 



Q(Q-2) + E2ii*2 



(18) 



Numerical data confirm that the isolated states are not 
affected by disorder up to a critical disorder strength, 
Wcr, that can be very large. Above W cr , the localization 
length of the isolated states also begins to decrease. This 



SUMMARY 

As an example of a non-trivial quantum network, we 
have studied the properties of an open circuit with M 
branches coupled to each other through one common 
point. The transmission properties of the open model 
were studied in relation to the structure of the energy 
spectrum and eigenfunctions of the parent closed system. 
The method of the effective non-Hermitian Hamiltonian 
allows one to derive the exact solution of the problem. 
Our main interest was in the energies and localization 
lengths of special eigenstates located outside the crystal 
energy band and strongly localized at the junction. It 
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tolica is also acknowledged. 



Figure 4: Critical value W C r in units where v = 1, as a function 
of Q = JZ a t 2 for M = 4. We choose for simplicity t a = t for 
a = 1 , . . . , M , so that Q = Mt 2 . Circles represent numerical data 
(full JV = 30, open N = 100), while the dashed line shows the 
prediction 1211 . 



was shown that, in the presence of coupling to the con- 
tinuum, these states typically become narrow resonances 
with a very large lifetime. This fact may be important for 
the fabrication of new kinds of electron nanostructures, 
waveguides, antennas, and lasing devices with a large 
quality factor. Another possible application of the M- 
branch circuits follows from the expression (|15j) . which 
points out the possibility of controlling and distribut- 
ing energy incoming via one branch into all others. We 
have also shown the negligible influence of disorder in the 
branches on these special eigenstates. In our considera- 
tion the continuum coupling 7 was taken as a constant 
parameter. In real arrangements, the circuit can be con- 
nected to transmitters or particle reservoirs. Then 7 can 
depend on the signal energy and density of states in the 
reservoirs. Along with that, the continuum coupling will 
acquire a real part (dispersive integral) that should be 
added to the intrinsic Hamiltonian in Eq. (10). We hope 
to consider this more complicated situation elsewhere. 
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